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CHAPTER I 


INTRODUCTION 

The primary objective of this report is to demonstrate the feasibility of using elec- 
trets, a new measuring device developed by the Space Sciences Laboratory of Marshall 
Space Flight Center (MSFC) to measure the chemical composition of rocket exhaust from 
space vehicles. The secondary objective is to evaluate the NASA/MSFC Multilayer Diffu- 
sion Model versus standard hydrogen chloride (HCl) detectors. It is noted there are over 
50 measuring techniques available to detect atmospheric contamination. The instrumenta- 
tion used operationally by the Booster Exhaust Study Test Team (BEST) which was com- 
posed of personnel from Langley Research Center (LaRC), Kennedy Space Center (KSC), 
and MSFC during the sampling of the ground cloud from the Titan IIIC launch of Viking 
I mission to Mars from Cape Kennedy, Florida, on August 20, 1975 and the instruments 
used for the measurement of gases in MSFC’s static test firings of the 6.4 percent scaled 
model of the Space Shuttle are presented. 

Several of these instruments have been used in previous attempts to obtain meas- 
urements of HCl concentrations downwind of actual launches of solid rockets, both at 
KSC and Vandenberg Air Force Base, California. However, due to the problems in trying 
to predict the path of the ground cloud and thereby prelocate the instruments, minimal 
data had been obtained. In such field tests, it soon becomes apparent that one must 
either deploy an extensive grid of measuring devices or have the instruments highly mobile 
and capable of following the cloud. It has been shown there is a general complication 
with all these types of instruments, i.e., a lack of multipollutant capabilities or a need for 
the development of a method of high sensitivity, simplicity, and speed of measurement, 
and for multipollutant capabilities. The electrets can overcome this complication. The 
experiments with the electrets have shown that they have unique properties that are use- 
ful in practical exhaust studies. They may be deployed over a large area and used to 
augment the more sophisticated detectors. 

The research is reported in the following manner. Chapter II discusses the back- 
ground, the preparation, and the use of electrets to measure effluents from Solid Rocket 



Boosters (SRB’s). Chapter III discusses the standard instruments used to detect HCl dur- 
ing the static tests, chamber tests, and Viking I mission. Chapter IV presents the meteoro- 
logical input which is vital to the NASA/MSFC Multilayer Diffusion Model to predict the 
HCl concentration. Chapter V discusses the NASA/MSFC Multilayer Diffusion Model. 
Chapters II, III, IV, and V present the necessary background information to evaluate real- 
istically the electrets. Chapter VI evaluates electrets by reporting the results of the static 
tests, chamber tests, and Viking I’s mission. Comparisons are made with the diffusion 
model and standard instruments. Finally, Chapter VII presents the conclusions of the 
investigation. 
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CHAPTER II 


EXPERIMENTAL ELECTRETS 

Electrets have been known since the latter part of the 19th century. A substantial 
amount of work has been done on electrets, and a few practical applications, particularly 
in the area of communications, have been made. In general, however, only qualitative 
descriptions of the various processes are known. 

The particles in tests that have been conducted have either a positive or negative 
charge and are attracted to the electret’s surface. These collected particles or ions on the 
electrets are then analyzed by taking the mass spectra of the sample. 

Electrets have been made in a variety of ways, but the method used most often 
has been to make “thermoelectrets.” A thermoelectret is made by placing a dielectric 
material between two electrodes, heating the dielectric material to near its softening 
point, and applying a high direct current field to the material between electrodes. The 
field is maintained while the dielectric is cooled to room temperature, and the electric 
field is removed. If the treated dielectric exhibits electric polarization, it is called an elec- 
tret [ 1 ] . 

It has been found that thermoelectrets of polymers with stable surface charges on 
either side are suitable devices for attracting charged particles and ions to their surface. 
Therefore, it was decided to use these electrets for effectively measuring charged gases, 
vapors, or particles in the atmosphere. The thermoelectret characteristics of polytetra- 
fluorethylene (Teflon) from experimental studies were analyzed by collecting vapors of 
xylene and acetone under different experimental conditions. Thermoelectrets were pre- 
pared by keeping a Teflon foil of 0.0125 cm thickness between two aluminum electrodes. 
Further details are in Reference 2. Mass spectra of the unexposed electret and the 
exposed electret were taken by using the Varian M-66 cycloidal mass spectrometer. The 
experiment results indicated that even small quantities of the gases give a characteristic 
mass spectra and the peak height of the mass spectra can be considered as a measure of 
the quantity of the material collected over the surface of the electret. 
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Electrets are being used to measure the effluents from the propellants in SRB’s. 
The most widely used propellant in SRB’s is ammonium perchlorate as the oxidizer with 
powdered aluminum filler that acts in part as a fuel and partially as a stabilizer to control 
the burning rate. The exhaust products from this type of fuel contain HCl, aluminum 
oxide (AI 2 O 3 ), and water (H 2 O). To assess the impact of these products in the atmos- 
phere, it is necessary to know not only their quantity but also their distribution in the 
ground cloud that develops at the launch site after a rocket firing. 

Electrets of polymers were used successfully for the first time in collecting rocket 
exhaust effluent measurements. The exposed electrets were analyzed by scanning electron 
microscope (SEM) and X-ray microprobe analyses. It was possible by these studies to 
identify the various effluents coming from the rockets at the time of firing. Direct com- 
parisons of values with electrets were made with the results obtained from the chemilumi- 
nescent, bubbler, coulometer, and millipore filters for samples kept under identical condi- 
tions. This aided in calibrating the intensity of the X-ray spectra with respect to the 
concentration of the products collected on the electret surfaces. The results of these com- 
parisons show a high degree of correlation. The entire spectra of the pollutants coming 
from the rocket may be obtained from the electrets. Other methods showed only the 
concentration of HCl. However, more extensive investigations must be performed to 
obtain quantitative results with electrets. This report discusses the use of electrets to 
measure the effluents. 

A typical computer output of the scanning electron microscope (SEM) and X-ray 
microprobe analysis [2] of the rocket exhaust effluents is shown in Figure 1 and listed in 
Table 1 . Figure 1 illustrates the X-ray energy spectra obtained from the electret during 
the Viking I launch (Test VI 1). Table 1 lists the peak count of X-ray energy for each 
element. The surface topology of the electret for this test is shown in Figure 2. Because 
this method has the added advantage of obtaining the entire spectrum of pollutants from 
a rocket exhaust, it is a powerful tool in studying contaminants. 
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i X-RAY ENERGY {KeV) 


Figure 1. X-ray energy spectra from Test VI 1, 
Viking I mission. 


Table 1. Data Output from Analyzer, Test VI 1, 
Viking I Mission 


ay Energy 
(KeV) 

Counts 
(2000 s) 

Symbol 

0.707 

13 632 

F 

1.497 

121 644 

A1 

1.749 

90 720 

Si 

2,105 

4 647 

Au 

2.306 

18 156 

S 

2.626 

24 992 

Cl 

3.317 

16 415 

K 

3.660 

3 146 

Ca 

5.409 

3 094 

Cr a 

5.880 

79 

Cr|3 

6.388 

1 

10 469 

Fe a 

7.035 

1 233 

Fep 

8.010 

24 390 

Cu a 



Figure 2. Surface topology from Test VII, 
Viking I mission. 
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CHAPTER III 


INSTRUMENTS 

The following instruments were used to detect and measure the concentration of 
HCl in these tests: 

1. Chemiluminescent HCl detector (Geomet) 

2. Bubblers 

3 . Millipore filter 

4. Coulometer. 

A brief discussion of each of the instruments and its evaluation follows. 


A. Chemiluminescent HCl Detector (Geomet) 


This detector was developed under contract for NASA and is now available com- 
mercially. The HCl detector is designed to continuously monitor HCl air concentrations 
from 0.1 to 50 ppm. In this concentration range, the response time to 90 percent signal 
and 10 percent recovery is a few seconds. Since these terms will be used again, an expla- 
nation is presented as follows [3] : 

Response time has been defined as that time required for the output of the 
detector to reach 90 percent of its final (equilibrated) output upon exposure 
to a fixed pollutant concentration. Time zero for a response time measure- 
ment is that time at which the pollutant concentration is introduced into the 
inlet of the detector. Recovery time is defined as that time required for the 
output of the detector sampling a fixed concentration to decay to 10 per- 
cent of the original value after removing source of the pollutant gas stream 
from the inlet of the detector. Time zero for a recovery time measurement is 
that time at which the source of pollutant is removed from the inlet of the 
detector. 

HCl is detected by chemiluminescence accompanying the reaction of luminol. The 
inlet column for the sampling stream is coated with NaBr and NaBrOa . HCl reacts with 
this coating to produce hypochlorite, hypobromite, and possibly bromine, all of which 
react with luminol to produce chemiluminescence which is proportional to the concentra- 
tion of HCl in the incoming stream. With most other instruments sufficient HCl, espe- 
cially at low concentrations, is removed by the adsorption or reaction at the inlet to 
cause a lag in response of the instrument. However, in this chemiluminescent instrument. 
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this adsorption or reaction produces species that are essential to the detection process. 
The response of this instrument is rapid in comparison to other HCl detecting instru- 
ments. Laboratory and field evaluations of the chemiluminescent HCl detector have 
shown the unit to be a rehable HCl detector [3] . Laboratory results have shown detector 
hmits of 0.05 ppm and an absolute accuracy of at least ±5 percent. Field results have 
shown similar detection limits and stable operation at ambient temperatures of 40° C. 
Future laboratory evaluations will be in the area of interference studies and detection of 
HCl (gas) and HCl (aerosol) mixtures. Based on the results in Reference 3, the following 
conclusions apply; 

1. The detection limit for the HCl detector is less than 0.05 ppm. 

2. In the range of 0.05 to 50 ppm, the accuracy of the HCl detector is 5 

percent. 

3. Response time for the HCl detector ranges from less than 1 s at 50 ppm to 
approximately 20 s at 0.05 ppm. 

B. Bubblers 

Basically, the bubbler is a dosage type instrument. It sucks air at a rate of 3000 
cm^/min into 20 cm^ of distilled water. It uses coulometric detection for the detection 
of chloride. 

There is a detection limit of 3 nanograms (ng)/20 pliters of water. Also, the 
detection limit is a function of how pure one can get the water. Basically, the detection 
limit is 80 to 100 ppm-s, which is marginal. As an example, if the concentration is 0.8 
ppm for 10 s, no measurement is obtained. However, if a 0.8 ppm measurement for 100 
s is made, a measurement may be obtained. The accuracy is 10 percent of the reading 
above 200 ppm-s; below this the accuracy increases [4] . 

C. Millipore Filter 

A millipore filter consists of a plastic membrane backed with a cellulose pad. The 
membrane and pad are housed in a disposable plastic container. The accuracy of the 
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determination of an HCl concentration using this device is defined by the various meas- 
urements that must be made, i.e., time of exposure, flow rate through filter, quantity of 
water used in the analysis, and the hydrogen ion concentration. Assuming reasonable 
uncertainties for each of these measurements yields a net uncertainty of ±15 percent in 
the HCl concentration. 

The pore size of the millipore filter membrane is 0.4 /rm. The flow rate to sample 
the exhaust gas is 5 liters/min. The membrane is 3.28 cm (1.5 in.) in diameter. The 
membrane collects the aluminum oxide particles. The HCl is absorbed by the membrane 
and the cellulose backing pad. The dosage absorbed by a specific filter is determined by 
removing the membrane and pad and macerating them in 20 cm^ of distilled water. The 
resulting pH change in the water is then measured. Thus, the hydrogen ion concentration 
can be determined, and the quantity of HCl follows from this determination [4] . 

D. Coulometer 

The coulometer used in these tests was supplied and operated by personnel from 
the U.S. Air Force School of Aerospace Medicine, Brooks Air Force Base, Texas. It is a 
standard laboratory microcoulometric titrating system- which has been modified to 
accept a continuous sample of gas bubbled through the electrolyte. The cell is composed 
of two pairs of electrodes immersed in a 70 percent (vol./vol.) aqueous glacial acetic acid 
solution made approximately 1 X 10''^ molar in silver ions. One pair of electrodes is 
used to sense the concentration of silver ions and the other pair is used to generate 
silver ions. When a sample of HCl is introduced into the cell, it reacts with the silver ions 
to produce AgCl as a precipitate. As the concentration of the silver ions in the electrolyte 
varies, a change in the output voltage across the sensing electrodes is developed. This volt- 
age is fed to an amplifier and after being amplified and properly phased is applied to the 
generating electrodes to replenish the silver ions. 

The potential drop across a precision resistor in series with the generating elec- 
trodes is monitored by a recorder. Thus, the total charge needed to regenerate the silver 
ions appears as a peak on the recorder. The area under this peak is related by Faraday’s 
Law to the CF ion concentration. The threshold detection limit for batch injection of 
HCl is approximately 3 nanograms (ng). The total quantity is calculated from 
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- 36-45 gm/mole ^ /a\ 

96 501 coulomb \^/ 

where 

W = weight of HCl nanograms (ng) 

A = coulogram peak area millivolts-seconds (mV-s) 

R = precision resistor ohms (f2) 

When the instrument is sampling a steady-state mixture of HCl and air, the result- 
ing constant current supplied to the generating electrodes can be directly related to the 
HCl concentration. The range of HCl concentrations which can be measured by the 
instrument can be adjusted by changing the rate at which the sampled gases are drawn 
through the electrolyte. However, there are limits, imposed by the excessive turbulence of 
the electrolyte caused by high flow rates and poor mixing at low sampling rates, which 
set the optimum instrument range from 0. 1 to 20 ppm of HCl. A typical sampling rate is 
0.1 liter/min. 

The coulometer has been employed as a standard to check HCl mixtures used as 
calibration gases for other instruments. In this role the coulometer samples the mixture 
for a specific time; from the determination of the weight of the total chloride ion sensed, 
the flow rate of the gas through the electrolyte, and the time interval, the concentration 
of the HCl is calculated. Using a chain rule analysis with estimates on the precision of 
each measured quantity indicates an expected uncertainty of ±6 percent when the cou- 
lometer is sampling at its upper limit of 17 ppm [4] . 
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CHAPTER IV 


METEOROLOGICAL DATA 

Meteorological data are used as weather input into the NASA/MSFC Multilayer 
Diffusion Model to compute the HCl isopleth concentration versus downwind distance 
and HCl centerline concentration as given in Chapter VI of this report. 

The atmospheric soundings for the Space Shuttle Vehicle’s (SSV’s) static test 
firings were made by the Atmospheric Research Test Facility of MSEC with the following 
equipment: 

a. Radiosonde — A balloon-borne instrument for the simultaneous measurement 
and transmission of meteorological data. The instrument consists of transducers for the 
measurement of pressure, temperature, and humidity. 

b. Rawin - A method of upper-air observation consisting of an evaluation of the 
wind speed and direction, temperature, and relative humidity aloft by means of a 
balloon-borne radiosonde tracked by a radar or by a ground meteorological detector. 
Height data pertaining to levels aloft are computed from the radiosonde data, while wind 
data are derived by. trigonometric computation. 

c. AF Type Radiosonde — At the KSC meteorological station, the rawin flights 
were made with an AF type radiosonde using the GMD-4 rather than the National 
Weather Service (NWS)/National Oceanic and Atmospheric Administration (NOAA) radio- 
sonde system to obtain weather data during Viking I’s launch. The temperature and 
humidity sensor data are transmitted ten times per minute in the AF sonde by a clock- 
actuated switch rather than the aneroid barometer switch used in the NOAA radiosonde. 
Both systems measure azimuth and elevation with the GMD. A transponder in the AF 
sonde is used to obtain the slant range to the radiosonde, enabling the calculation of alti- 
tude. The pressure is then calculated according to the hydrostatic equation. The equa- 
tions used in the computer program to calculate various thermodynamic quantities from 
the basic altitude, temperature, and relative humidity data are given in Appendix A. 
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A. Static Tests’ Weather Data for Diffusion Model 


Tables 2 through 19 present a collection of atmospheric soundings of the environ- 
ment at MSFC’s Test FacUity 116, Huntsville, Alabama, during the SSV 6.4 percent 
model static test firings. 

B. Viking Fs Weather Data for Diffusion Model 

Table 20 is the atmospheric soundings for the Viking I launch. The Viking I 
mission to Mars was a Titan HIE launch from Pad 41 at KSC at 1722 EDT on August 
20, 1975. 

The test day occurrences of the 18 static test firings were as follows: two in 
August, two in November, one in December, nine in January, three in February, and one 
in March. The nine tests in January were particularly favorable for model testing of HCl 
downwind concentrations because of the extreme wind shears, temperature inversions, 
and cold temperatures. 


Table 2. GMD Sounding for Test Firing No. 1, August 16, 1974, 
1440 CSX MSFC Test Facility 116^ 


Height 

(m) 

Wind Speed 
(m/s) 

Wind Direction 
(degree) 

Temperature 

(°C) 

Pressure 

(mb) 

0.0 

2.6 

240.0 

29.3 

999.0 

107.0 

3.6 

229.8 

26.9 

987.0 

197.0 

4.6 

219.6 

26.9 

977.0 

306.0 

4.0 

207.7 

25.9 

965.0 

434.0 

4.4 

210.5 

24.5 

951.0 

573.0 

4.2 

214.6 

22.4 

936.0 

695.0 

4.1 

209.6 

21.3 

923.0 

837.0 

4.2 

208.7 

20.0 

908.0 

971.0 

3.9 

208.9 

19.4 

894.0 

1087.0 

4.1 

209.6 

18.5 

882.0 


a. Standard deviation of the surface wind angle — 4.0°; surface air density — 
1156.0 gm/m^ 



Table 3. GMD Sounding for Static Test Firing No. 2, August 30, 1974, 
1302 CST, MSFC Test Facility 116^ 


Height 

(m) 

Wind Speed 
(m/s) 

Wind Direction 
(degree) 

Temperature 

(°C) 

Pressure 

(mb) 

0.0 

3.6 

250.0 

25.3 

997.9 

123.0 

4.6 

246.5 

24.4 

984.0 

221.0 

5.7 

243.0 

23.4 

973.0 

357.0 

5.9 

238.4 

22.1 

958.0 

485.0 

6.6 

236.0 

21.0 

944.0 

679.0 

8.3 

234.0 

19.2 

923.0 

897.0 

8.9 

230.6 

18.2. 

900.0 

1090.0 

9.4 

222.4 

17.8 

880.0 


a. Standard deviation of the surface wind angle — 4.0°; surface air density — 
1218.9 gm/m^ 


Table 4. GMD Sounding for Static Test Firing No. 3, November 19, 1974, 
2131 CST, MSFC Test Facility 116^ 


Height 

(m) 

Wind Speed 
(m/s) 

Wind Direction 
(degree) 

Temperature 

(°C) 

Pressure 

(mb) 

0.0 

0.5 

100.0 

13.6 

989.0 

170.0 

3.9 

153.6 

14.3 

979.0 

284.0 

7.3 

207.2 

14.9 

965.0 

408.0 

7.6 

219.2 

15.1 

951.0 

542.0 

8.3 

278.3 

13.5 

936.0 

687.0 

9.6 

236.2 

12.1 

920.0 

834.0 

10.1 

243.2 

10.6 

904.0 

973.0 

10.6 

244.7 

9.6 

889.0 

1105.0 

11.6 

250.3 

8.3 

875.0 


a. Standard deviation of the surface wind angle — 4.0°, surface air density 
1215.95 gm/m^ . 
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Table 5. GMD Sounding for Static Test Firing No. 4, November 22, 1974, 
1830 CST, MSFC Test Facility 1 16^ 


Height 

Wind Speed 

Wind Direction 

Temperature 

Pressure 

(m) 

(m/s) 

(degree) 

(°C) 

(mb) 

0.0 

0.5 

90.0 

7.1 

1000.0 

101.0 

4.0 

127.6 

16.0 

988.0 

248.0 

7.6 

165.2 

15.1 

971.0 

407.0 

4.6 

157.0 

14.4 

953.0 

523.0 

2.0 

162.9 

13.6 

940.0 

649.0 

1.8 

196.2 

12.6 

926.0 

745.0 

1.4 

243.4 

12.2 

910.0 

925.0 

2.6 

276.2 

12.6 

896.0 

1057.0 

3.4 

270.9 

12.8 

882.0 


a. Standard deviation of the surface wind angle — 4.0°; surface air density — 
1244.4 gm/m^ . 


Table 6. GMD Sounding for Static Test Firing No. 5, December 12, 1974, 
1350 CST, MSFC Test Facility 1 16^ 


Height 

Wind Speed 

Wind Direction 

Temperature 

Pressure 

(m) 

(m/s) 

(degree) 

(°C) 

(mb) 

0.0 

4.1 

190.0 

17.1 

997.9 

136.0 

5.2 

187.1 

14.4 

982.0 

265.0 

6.3 

184.2 

13.2 

967.0 

397.0 

6.9 

191.2 

12.0 

952.0 

511.0 

7.5 

198.6 

11.1 

939.0 

627.0 

8.5 

199.4 

9.6 

926.0 

736.0 

9.4 

204.0 

8.8 

914.0 

845.0 

10.2 

208.2 

7.9 

902.0 

955.0 

10.7 

211.0 

7.7 

890.0 

1077.0 

10.0 

213.3 

7.9 

877.0 


a. Standard deviation of the surface wind angle — 4.0°; surface air density — 
1197.7 gm/m^. 
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Table 7. GMD Sounding for Static Test Firing No. 6, January 7, 1975, 
1759 CST, MSFC Test Facility 1 16 ^ 


Height 

(m) 

Wind Speed 
(m/s) 

Wind Direction 
(degree) 

Temperature 

(°C) 

Pressure 

(mb) 

0.0 

2.6 

130.0 

12.7 

997.0 

152.0 

6.9 

144.4 

12.4 

979.0 

291.0 

11.3 

158.9 

12.0 

963.0 

404.0 

12.2 

165.8 

11.2 

950.0 

528.0 

12.5 

172.0 

10.8 

936.0 

645.0 

13.5 

177.1 

10.1 

923.0 

772.0 

13.3 

182.8 

9.7 

909.0 

919.0 

12.8 

189.8 

9.1 

893.0 

1069.0 

12.4 

197.5 

8.4 

877.0 


a. Standard deviation of the surface wind angle — 4.0°; surface air density — 
1218.9 gm/m^ 


Table 8. GMD Sounding for Static Test Firing No. 7, January 9, 1975, 
1758 CST, MSFC Test Facility 1 16^ 


Height 

(m) 

Wind Speed 
(m/s) 

Wind Direction 
(degree) 

Temperature 

(°C) 

Pressure 

(mb) 

0.0 

5.2 

110.0 

14.3 

991.5 

141.0 

9.3 

129.0 

13.9 

975.0 

290.0 

13.5 

148.1 

13.5 

958.0 

450.0 

13.8 

162.8 

13.2 

940.0 

604.0 

13.2 

181.4 

14.1 

923.0 

733.0 

12.8 

198.3 

15.1 

909.0 

846.0 

13.3 

208.4 

14.7 

897.0 

979.0 

14.3 

212.4 

14.1 

883.0 

1114.0 

15.3 

215.2 

13.2 

869.0 


a. Standard deviation of the surface wind angle — 4.0° ; surface air density 
1199.12 gm/m^ . 
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Table 9. GMD Sounding for Static Test Firing No. 8, January 1 1, 1975, 
1435 CST, MSFC Test Facility 1 16^ 


Height 

Wind Speed 

Wind Direction 

Temperature 

Pressure 

(m) 

(m/s) 

(degree) 

CC) 

(mb) 

0.0 

2.6 

230.0 

12.5 

996.4 

198.0 

3.0 

255.0 

10.3 

973.0 

396.0 

3.4 

280.1 

8.4 

950.0 

571.0 

3.4 

271.7 

6.4 

930.0 

749.0 

4.2 

264.5 

4.3 

910.0 

939.0 . 

4.9 

257.3 

2.8 

889.0 

1113.0 

5.4 

246.6 

0.9 

870.0 


a. Standard deviation of the surface wind angle — 4.0°; surface air density — 
1212.64 gm/m^ . 


Table 10. GMD Sounding for Static Test Firing No. 9, January 15, 1975, 
1550 CST, MSFC Test Facility 1 16a 


Height 

(m) 

Wind Speed 
(m/s) 

Wind Direction 
(degree) 

Temperature 

(°C) 

Pressure 

(mb) 

0.0 

1.0 

270.0 

10.3 

1001.6 

130.0 

2.0 

252.9 

8.9 

986.0 

273.0 

3.1 

238.8 

7.8 

969.0 

410.0 

4.7 

241.7 

6.2 

953.0 

540.0 

6.5 

242.6 

5.3 

938.0 

698.0 

7.7 

246.5 

4.5 

920.0 

841.0 

7.9 

261.0 

4.3 

904.0 

968.0 

8.1 

276.1 

4.1 

890.0 

1106.0 

8.7 

287.3 

2.8 

875.0 


a. Standard deviation of the surface wind angle — 4.0°; surface air density — 
1230.60 gm/m^ . 
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Table 11. GMD Sounding for Static Test Firing No. 10, January 16, 1975, 
1640 CST, MSFC Test Facility 1 16^ 


Height 

(m) 

Wind Speed 
(m/s) 

Wind Direction 
(degree) 

Temperature 

(°C) 

Pressure 

(mb) 

0.0 

1.6 

330.0 

7.2 

1006.4 

135.0 

4.7 

350.1 

6.8 

990.0 

260.0 

7.8 

10.3 

5.8 

975.0 

361.0 

7.5 

12.8 

4.6 

963.0 

472.0 

7.3 

11.7 

3.4 

950.0 

566.0 

7.2 

16.4 

2.2 

939.0 

678.0 

7.8 

16.7 

1.3 

926.0 

810.0 

7.7 

16.0 

0.4 

911.0 

943.0 

6.6 

11.8 

-0.8 

896.0 

1069.0 

5.0 

359.5 

-1.3 

892.0 


a. Standard deviation of the surface wind angle - 4.0°; surface air density — 
1249.88 gm/m^. 


Table 12. GMD Sounding for Static Test Firing No. 1 1, January 18, 1975, 
1450 CST, MSFC Test Facility 1 16^ 


Height 

(m) 

Wind Speed 
(m/s) 

Wind Direction 
(degree) 

Temperature 

(°C) 

Pressure 

(mb) 

0.0 

0.5 

350.0 

1 1.7 

996.3 

138.0 

4.4 

294.2 

1 1.9 

980.0 

275.0 

8.4 

238.5 

11.3 

964.0 

415.0 

10.7 

246.7 

10.6 

948.0 

557.0 

13.2 

248.8 

9.7 

932.0 

701.0 

16.4 

247.7 

8.9 

916.0 

847.0 

17.5 

252.8 

7.7 

900.0 

986.0 

17.7 

256.4 

7.5 

885.0 

1126.0 

18.2 

260.1 

6.6 

870.0 


a. Standard deviation of the surface wind angle — 4.0° ; surface air density — 
1245.46 gm/m^ . 
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Table 13. GMD Sounding for Static Test Firing No. 12, January 21, 1975, 
1 145 CST, MSFC Test Facility 1 16^ 


Height 

Wind Speed 

Wind Direction 

Temperature 

Pressure 

(m) 

(m/s) 

(degree) 

(°C) 

(mb) 

0.0 

3.1 

120.0 

7.5 

1004.0 

132.0 

4.7 

130.4 

7.3 

988.0 

266.0 

6.3 

140.9 

6.0 

972.0 

402.0 

7.6 

150.3 

4.8 

956.0 

539.0 

7.6 

155.7 

2.7 

940.0 

678.0 

6.9 

165.3 

1.4 

924.0 

801.0 

6.5 

177.2 

1.6 

910.0 

944.0 

6.7 

192.2 

3.2 

894.0 

1081.0 

7.3 

308.8 

3.6 

879.0 


a. Standard deviation of the surface wind angle — 4.0° ; surface air density — 
1249.87 gm/m^ 


Table 14. GMD Sounding for Static Test Firing No. 13, January 23, 1975, 
1958 CST, MSFC Test Facility 1 16^ 


Height 

Wind Speed 

Wind Direction 

Temperature 

Pressure 

(m) 

(m/s) 

(degree) 

(°C) 

(mb) 

0.0 

0.5 

40.0 

6.0 

1003.8 

139.0 

1.4 

80.7 

9.5 

987.0 

266.0 

2.3 

121.4 

8.7 

972.0 

403.0 

2.5 

125.0 

7.6 

956.0 

524.0 

3.0 

144.3 

7.3 

942.0 

639.0 

3.7 

163.6 

6..4 

929.0 

736.0 

5.1 

180.0 

8.0 

915.0 

900.0 

6.7 

198.0 

9.1 

900.0 

1039.0 

8.1 

213.7 

9.3 

885.0 


a. Standard deviation of the surface wind angle - 4.0°; surface air density — 
121 5.70 gm/m^. 
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Table 15. GMD Sounding for Static Test Firing No. 14, January 27, 1975, 
1822 CST, MSFC Test Facility 1 16^ 


Height 

(m) 

Wind Speed 
(m/s) 

Wind Direction 
(degree) 

Temperature 

(°C) 

Pressure 

(mb) 

0.0 

2.6 

160.0 

18.2 

996.6 

144.0 

7.1 

174.3 

19.5 

980.0 

259.0 

11.6 

188.7 

19.2 

967.0 

393.0 

14.5 

190.8 

18.6 

952.0 

529.0 

15.2 

196.9 

16.4 

937.0 

656.0 

16.8 

202.0 

15.0 

923.0 

795.0 

17.8 

209.2 

14.6 

908.0 

936.0 

18.1 

219.2 

14.2 

893.0 

1070.0 

19.3 

229.7 

14.8 

879.0 

a. Standard deviation of the surface wind angle — 4.0°; surface air density — 
1195.11 gm/m^. 

Table 16. GMD Sounding for Static Test Firing No. 15, February 8, 1975, 
1254 CST, MSFC Test Facility 1 16^ 

Height 

(m) 

Wind Speed 
(m/s) 

Wind Direction 
(degree) 

Temperature 

(°C) 

Pressure 

(mb) 

0.0 

2.1 

220.0 

7.8 

1000.2 

158.0 

2.1 

215.1 

3.7 

981.0 

325.0 

2.2 

210.3 

1.8 

961.0 

485.0 

4.0 

217.5 

0.3 

942.0 

657.0 

5.4 

231.0 

-1.4 

922.0 

823.0 

6.3 

236.0 

-1.5 

903.0 

938.0 

7.8 

241.2 

-0.4 

890.0 

1056.0 

8.3 

243.6 

-0.7 

877.0 


a. Standard deviation of the surface wind angle — 4.0°; surface air density 
1243.70 gm/m^ . 
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Table 17. GMD Sounding for Static Test Firing No, 16, February 12, 1975, 
1700 CST, MSFC Test Facility 1 16^ 


Height 

(m) 

Wind Speed 
(m/s) 

Wind Direction 
(degree) 

Temperature 

(°C) 

Pressure 

(mb) 

0.0 

2.6 

340.0 

5.8 

999.0 

106.0 

4.9 

357.0 

4.0 

986.0 

214.0 

7.2 

14.1 

2.5 

973.0 

323.0 

7.4 

17.4 

1.5 

960.0 

433.0 

7.3 

13.3 

0.6 

947.0 

543.0 

8.2 

19.0 

-0.3 

934.0 

673.0 

8.8 

22.5 

-1.3 

919.0 

795.0 

9.6 

30.4 

-2.2 

905.0 

929.0 

10.3 

39.3 

1.2 

890.0 

1048.0 

11.0 

44.7 

5.2 

877.0 

a. Standard deviation of the surface wind angle — 4.0°; surface air density — 
1251.30 gm/m^ 

Table 18. GMD Sounding for Static Test Firing No. 17, February 22, 
1528 CST, MSFC Test Facility 1 16^ 

1975, 

Height 

Wind Speed 

Wind Direction 

Temperature 

Pressure 

(m) 

(m/s) 

(degree) 

(°C) 

(mb) 

0.0 

7.2 

130.0 

13.8 

993.8 

151.0 

10.0 

134.4 

11.7 

976.0 

298.0 

12.9 

138.8 

10.2 

959.0 

447.0 

18.1 

154.0 

9.9 

942.0 

590.0 

20.8 

169.3 

11.4 

926.0 

700.0 

23.1 

178.9 

14.1 

914.0 

831.0 

26.3 

184.8 

14.7 

900.0 

973.0 

25.9 

191.1 

13.7 

885.0 

1097.0 

25.5 

195.9 

12.3 

872.0 


a. Standard deviation of the surface wind angle — 4.0°; surface air density — 
1210.02 gm/m^. 
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Table 19. GMD Sounding for Static Test Firing No. 18, March 25, 1975, 
1700 CST, MSFC Test Facility 1 16^ 


Height 

(m) 

Wind Speed 
(m/s) 

Wind Direction 
(degree) 

Temperature 

(°C) 

Pressure 

(mb) 

0.0 

6.7 

300.0 

12.9 

998.0 

177.0 

7.5 

302.7 

11.6 

977.0 

350.0 

8.8 

305.5 

9.7 

957.0 

524.0 

9.4 

298.3 

7.9 

937.0 

702.0 

10.2 

294.1 

5.9 

917.0 

882.0 

11.7 

286.6 

3.8 

897.0 

1037.0 

13.6 

283.7 

2.0 

880.0 


a. Standard deviation of the surface wind angle — 4.0°; surface air density — 
1215.8 gm/m^. 


Table 20. Titan T-0 Sounding, August 20, 1975, 
1722 EDT(2122Z)^ 


Height 

(m) 

Wind Speed 
(m/s) 

Wind Direction 
(degree) 

Temperature 

(°C) 

Pressure 

(mb) 

5.0 

3.0 

1 10.0 

28.7 

1018.3 

166.0 

4.0 

95.0 

26.8 

1000.0 

220.0 

4.0 

91.0 

26.2 

994.0 

305.0 

4.0 

89.0 

25.5 

984.5 

610.0 

2.0 

83.0 

24.5 

951.1 

619.0 

2.0 

82.0 

24.5 

950.0 

638.0 

2.0 

82.0 

24.5 

948.0 

814,0 

1.0 

87.0 

21.8 

918.5 

1078.0 

0.1 

77.0 

20.7 

902.0 

1090.0 

0.1 

80.0 

20.6 

900.0 

1219.0 

0,1 

109.0 

20.0 

886.9 

1524.0 

1.0 

169.0 

17.7 

856.1 

1583.0 

1.0 

166.0 

17.4 

850.0 

1829.0 

1.0 

162.0 

16.8 

826.1 

1987.0 

3.0 

167.0 

15.4 

81 1.0 

2100.0 

4.0 

170.0 

15.5 

800.0 

2134.0 

4.0 

170.0 

15.6 

797.1 

221 1.0 

5.0 

171.0 

15.6 

790.0 

2438.0 

4.0 

159.0 

1 1.7 

769.0 


a. Surface air density - 1 162.670 g/m’ ; height of surface mixing layer 
1090 m. 
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CHAPTER V 


NASA/MSFC MULTILAYER DIFFUSION 
MODELS PROGRAM 

The NASA/MSFC Multilayer Diffusion Models Program is designed to calculate 
the following quantities downwind from normal and abnormal launches of rocket vehi- 
cles: 

a. Concentration and dosage patterns 

b. Time interval (minutes) — concentration 

c. Average cloud concentration 

d. Time of cloud passage 

e. Ground-level deposition patterns due to gravitational settling or precipitation 
scavenging. 

Programs options include the calculation of concentration, dosage, and time-mean con- 
centration patterns with partial reflection of material at the surface, with time-dependent 
exponential decay, and/or with depletion due to precipitation scavenging. Also, the pro- 
gram is capable of calculating ground-level deposition due to gravitational settling with 
partial reflection at the surface. Provision is also made (in Model 4) to account for 
changes in meteorological structure along the cloud trajectory. Model 4 also may be used 
to determine concentration and dosage fields in the surface mixing layer downwind from 
a source in which the source strength varies with height in the layer [5] . 

Program output options include: 

a. Printing of all data inputs 

b. Printing of the results of all model calculations 

c. Plotting of maximum centerline concentration, dosage, time-mean concentra- 
tion, and deposition versus distance along the cloud trajectory 

d. Plotting of concentration, dosage, time-mean concentration, and deposition 
isopleths. 

A simplified block diagram illustrating major features of the program is shown in 
Figure 3. 
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With the increasing knowledge of the environmental impact of air pollution, 
special attention is now being placed on the emission into the atmosphere of aerospace 
vehicle exhaust effluents and by-products. Limited concern was placed on this problem in 
the past, especially because vehicles using solid rocket motors were small in size and few 
were tested and launched. National, state, and local air pollution laws are also becoming 
more stringent to prevent environmental damage which places additional restrictions on 
any organization or private citizen in regard to pollution of many types. 

To determine the ground level concentrations, the procedure is to apply the 
proper atmospheric diffusion models to calculate downwind concentrations and dosages 
from various engine and SRB exhaust by-products. A major effort is being made to 
gather detailed data on the chemical reactions that take place between the exhaust efflu- 
ents and the atmosphere. While little is currently known about this problem, research is 
underway throughout the aerospace research community to determine initial and long 
term source characteristics. 
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This chapter includes statements on the basic diffusion estimation formulas. This 
is a summary of the salient facets of the effluent material transport problem which is 
found in referenced hterature. Other than normal exhaust releases and abnormal releases, 
brief statements will follow on leaks and inadvertent spills. Cloud rise formulas for use in 
source identification are included and meteorological inputs are covered. 

A. Definitions 

Concentration is the mass of a pollutant per unit volume at a point in space and 
is referenced to the ambient atmosphere (units: parts per milhon, milligrams per cubic 
meter, etc.). 

Dosage is the time-integrated concentration at a point in space and has the units 
of concentration multiplied by the unit of time. 

The generalized concentration and dosage models describe the behavior of the 
cloud of toxic material after the cloud establishes equilibrium in the mixing layer. This 
equilibrium point is known as the cloud source and serves as the origin for the Cartesian 
coordinate system such that the x-axis is in the mean azimuth wind direction, y is the 
crosswind or lateral direction, and z is the vertical height above the ground. (The location 
of the source cloud is addressed later.) It is also assumed that this is an expanding volu- 
metric cloud about a moving reference point in a homogeneous fluid. For diagnostic and 
interpretation flexibility, these models are formatted in a modular form [5,6]. 

The generalized concentration model for a nearly instantaneous source is 
expressed as the product of seven modular terms: 

Concentration = (Peak Concentration’ Term } X {Alongwind Term} 

X (Lateral Term } X (Vertical Term} X (Depletion Term} 

X (Scavenging Term} X (Surface Absorption Term} 

whereas, the generalized dosage model for a nearly instantaneous source is defined by the 
product of six modular terms: 
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Dosage = {Peak Dosage Terms } X {Lateral Term} 

X {Vertical Term } X {Depletion Term} 

X {Scavenging Term } X { Surface Absorption Term } 

Thus, the mathematical descripfion for the concentration and dosage models per- 
mitsflexibility in application to various sources and for changing atmospheric parameters 
while always maintaining a rigorous mass balance. 

Two obvious differences exist. First, the peak concentration term refers to the 
concentration at the point x, y = O, z = H and is defined by the expression 

Q 

Point Peak Concentration = 3 , (5.1) 

(27t) ‘^x^y'^z 

where Q is the source strength and oj is the standard deviation of the concentration dis- 
tribution in the ith direction; whereas, the peak dosage term is given by 


Point Peak Dosage = 


Q 

27Tuaya2 


(5.2) 


where u is the mean wind speed. The second difference between these models is that the 
concentration contains a modular alongwind term to account for downstream temporal 
effects not considered in the dosage model. The alongwind term affords an exponential 
decay in concentration as a function of cloud transit time, concentration distribution, 
and the mean wind speed. 

The lateral term, common to both models, is another exponential decay term and 
is a function of the Gaussian spreading rate and the distance laterally from the mean 
wind azimuth. The vertical term, again common to both models, is a rather complex 
decay function since it contains a multiple reflection term for the point source which 
stops the vertical cloud development at the top of mixing layer and eventually changes 
the form of the vertical concentration distribution from Gaussian to rectangular. The last 
remaining three terms represent the options. The deposition term accounts for gravita- 
tional settling. The scavenging term accounts for precipitation of effluents by rain falling 


25 



through the exhaust cloud. The surface absorbtion term accounts for the fraction of 
material absorbed at a surface. 

The ground cloud is formed during the first several minutes and contains all of 
the exhaust by-products formed by the rocket engines from the time of engine ignition 
until the vehicle passes through the stabilization height, “height of maximum buoyancy 
rise of the hot exhaust products,” of the ground cloud. If the stabilization height is such 
that some of the ground cloud is contained in a thermally stable layer, only a fraction of 
the total amount of exhaust products in the ground cloud is available for mixing to the 
ground surfaces. 

The exhaust trail plume is the plume of stabilized exhaust products formed by 
rocket engine emissions occurring above the stabilization height of the ground cloud. 

B. Description of Models 

The normal launch environment will usually involve an atmospheric structure 
comprised of several horizontal meteorological layers with distinctive wind velocity, tem- 
perature, and humidity regimes between the surface and 5 km altitude. Large horizontal 
spatial variation in these meteorological parameters may also occur in the surface layer as 
a consequence of changes in terrain or land-water interfaces, which is accounted for by 
the absorbtion coefficient. The general diffusion model for concentration [equation 
(5.1)] and the dosage [equation (5.2)] assumes an expanding volume about a moving 
point of reference in a homogeneous environment. 

To overcome the obvious shortcomings of the general diffusion model but to stay 
within the accepted bounds of classical fluid mechanics [6] , a multiple layer concept is 
introduced to cope with the vertical and horizontal atmospheric gradients. Here, the gen- 
eral diffusion model is applied to individual horizontal layers in which the meteorological 
structure is reasonably homogeneous and independent of the neighboring layers. These 
layers have boundaries which are placed at points of major discontinuities in the vertical 
profiles of wind velocity, temperature, and humidity. Since the NASA/MSFC Multilayer 
Diffusion Model [5] is an empirical model and has imposed the general restriction of 
layer independence (no flux of particles or gases entering or leaving an individual layer), 
special provision must be made for spatial changes in the horizontal meteorology and for 
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gravitational settling or precipitation scavenging. In addition, the type of source within a 
layer must be considered, that is, whether there is a ground cloud source or a plume 
cloud source. 

The NASA/MSFC Multilayer Diffusion Model has six submodels to deal with the 
stages of the development of the exhaust cloud and the complex potentially varying 
meteorological conditions. These submodels can be used alone to describe all the environ- 
mental layers or in combinations where variations in layer meteorology require different 
modeling. For the introductory overview, however, these combinations will not be con- 
sidered. The primary output of all submodels is a mapping of the regimes of the con- 
centration and dosage isopleths. 

Model 1 is the basic model for the dispersive description of exhaust material from 
rocket plume. In this model it is assumed that the source extends vertically through the 
entire layer with a uniform distribution of the concentration of exhaust material, 
whereas, the horizontal distribution of the material being dispersed along the layer (x and 
y directions) has a Gaussian distribution. In addition, it is assumed that there is turbulent 
mixing. 

An analogous model would be the wave generated by dropping a rock in a river, 
where the wave disperses across the surface of the moving river. The significance of the 
supposition of turbulent mixing is that this mixing action disperses the effluent material 
across the layer similar to the way the wave is dispersed across the surface of the water. 
This model is an effective description of the plume cloud where the action of the vehicle 
passing through a layer leaves a cylindrical cloud of exhaust effluents. 

Model 2 is the same as Model 1 except it is assumed that there is no turbulent 
mixing. This implies that the exhaust material meanders along the layer without dis- 
persing, very much like a small oil puddle moves on the surface of a river. While Model 2 
is not generally used, movies of rocket firings clearly show that under some special 
meteorological conditions (temperature inversion) this model is required. While the multi- 
layer diffusion model is general in applicability, it is specific in meteorological parameters 
and launch description. 

From the standpoint of environmental impact, the description of the fields of the 
ground deposition of materials from the ground cloud is of primary significance — this 
description is afforded by Model 3. Generally, this model is employed in the surface 
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layer, but can be employed in any layer where the source does not extend through the 
entire layer. In this model, a Gaussian distribution is assumed along all three axes with 
turbulent mixing occurring. 

Model 3 is similar to Model 1 except that, rather than the plane two-dimensional 
dispersion of Model 1 , there is a three-dimensional dispersion of the exhaust cloud as the 
cloud is transported downstream. When the cloud reaches the top of the mixing layer, 
the distribution of material is reflected back into the expanding vertical distribution. 

Thus, after a certain time period Model 3 is identical with Model 1 . A clear understand- 
ing can be obtained if the formulation for this model is examined. 

The first three models can be summarized as describing initial transport of the 
effluents after the cloud reaches equilibrium. While the equations given for Model 3 are 
in the general form for any Kth layer, it should be noted that K = 1 (the surface layer) is 
used in most applications of this model. If after launch the rocket explodes in a layer, 
this can be studied with this model by setting K equal to that layer number. 

The remaining models are specialized models which afford a second order descrip- 
tion of the transport of the vehicle exhaust materials. These three models incorporate 
considerations for changes in meteorological conditions and particle effects. 

Model 4 updates the diffusion model with changes in meteorological conditions 
and structure which can occur as the toxic cloud propagates downstream. This model 
assumes that the vertical concentration of material has become uniform throughout each 
layer when a step-change in the meteorological conditions is introduced. This step-change 
results in the destruction of the original layer boundaries and the formation of new layer 
boundaries. The concentration fields which exist at this time are treated as new sources. 

In those new layers which now comprise more than one old layer, the old concentration 
is mapped as two independent concentration sources and then superimposed for the 
resulting concentration and dosage mappings. 

Model 5 accounts for precipitation scavenging. An example of where Model 5 
must be used is in solid rocket launches during the occurrence of rain, because the HCl 
will be scavenged by the rain. 

Model 6 describes the ground deposition due to gravitational settling of particles 
of droplets. Wind shears are incorporated in this model to account for the effect of the 
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settling velocity of the particulate matter. There are two forms for the source in this 
model, namely; 

1 . The source that extends vertically through the entire layer with a uniform 
distribution (this is the same source model as used with Models 1 and 2) 

2. A volume source in the Kth layer (this is the same source model as used with 
Model 3). 

Model 6 is very important in the analysis of the settling of AI2O3 particles 
released in solid rocket firings. 

In summary, the Multilayer Diffusion Model is composed of six submodels. 

Models 1 and 3 are designed to distinguish between the two sources of toxic cloud 
formation — the ground cloud during the initial launch phase (Model 3) and the plume 
cloud after the initial launch phase (Model 1). Model 2 was injected to account for a lack 
of turbulent mixing which can occur in the upper atmosphere. Model 4 is employed 
when a change in meteorological condition occurs during the downstream transport of 
the cloud. In the event of rain, the precipitation scavenging, both of gases and particles, 
can be accounted for in Model 5. The fallout of particulate matter on the ground is the 
domain of Model 6. These six submodels form the basic algorithms which are available to 
treat the diffusion problem. To model a specific launch of a vehicle, it is necessary to 
blend these algorithms together and adjust the model parameters to the specific meteoro- 
logical conditions of the launch, to the specific terrain (use of absorbtion coefficient) 
around the launch site, and to the specific vehicle being launched; this is the degree of 
complexity in the diffusion model. 

C. Cloud Rise Calculations 

The burning of rocket engines during normal launches results in the formation of 
a cloud of hot exhaust products which subsequently rises and entrains ambient air until 
an equilibrium with ambient conditions is achieved. Experience in predicting the buoyant 
rise [7] from normal launches of solid-fueled vehicles indicates the rise is best predicted 
using a cloud rise model for instantaneous sources. For solid-fueled vehicles, residence 
times near the pad are relatively short. 

Each of the models for cloud height is subdivided into two categories to account 
for the atmospheric temperature lapse rate. The model assumes that the atmosphere is 
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either quasi-adiabatic or stable. Here the quasi-adiabatic is where the adiabatic atmosphere 
is the limit, which means that the potential temperature difference (A0) is zero or less, 
where the potential temperature difference is given by 

Ad - cloud height ' ^surface • (5-3) 

If this potential temperature difference is positive, then the atmosphere is stable. Since 
in most cases of interest there will be an inversion layer present, the stable cloud rise 
formula is the normally utilized relation [8]. 

The maximum rise Zj^j for an instantaneous source is given by the expression in a 
stable atmosphere as 
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(5.4) 


whereas the maximum cloud rise downwind from an instantaneous source in a 
quasi-adiabatic atmosphere is given by 


^ml ~ 
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(5.5) 


In deriving equation (5.4), it is assumed that the initial upward momentum imparted to 
the exhaust gases by reflection from the ground surface and launch pad hardware is 
insignificant in comparison with the effect of thermal buoyancy. Based on limited experi- 
ence in predicting cloud rise from launches at Vandenberg Air Force Base, California, 
this assumption appears to be justified. The time required for the cloud to reach the 
stabilization height is given by the expression 

7T 




(5.6) 


In calculating z^^j from equations (5.4) and (5.5), the instantaneous heat released 
Qj is obtained from the relationship 


Ql - QptR (Zml } 


(5.7) 
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Inspection of the previously given equations reveals an interdependence between 
the calculated maximum cloud rise Zj^j, the height over which the potential temperature 
gradient 9<I>/9z is measured, and the value of used in obtaining Qj. Thus, the 

final value of maximum cloud rise must be found through iteration of either equation 
(5.4) or (5.5). The height over which 9<h/9z is measured and the time tR{ Zj^l are thus 
made consistent with the value of Zj^j calculated from the model. 


D. Cloud Dimensions and Vertical Distribution 
of Exhaust Products 


Source inputs required for the diffusion model calculations include the stabiliza- 
tion height of the exhaust cloud and cloud dimensions, as well as the vertical distribution 
of exhaust products in the stabilized cloud. The calculation of the stabilization height Zj^^ 
was described on page 29. The calculation of the dimensions of the stabilized cloud and 
the vertical distribution of exhaust products is described in the following paragraphs. 

The general formula used to calculate the radius of the stabilized cloud at height 
z is given by the expression 


(z) 


7Z 


; z 




-m 


y(2Zj^ - z) > 200 m; z > z 


m 


(5.8) 


where 


7 ~ (7i) " 0.64 

Note that for z > Zj^, the minimum radius of the stabilized cloud is set equal to 200 m. 

The cloud is assumed to be symmetrical about a vertical axis through the cloud 
centroid. The alongwind and crosswind source dimensions of the cloud in each of the 
layers are calculated under the following assumptions: 

1. The distribution of exhaust products within the cloud is Gaussian in the hori- 
zontal plane. 

2. The concentration of exhaust products at a lateral distance of one radius from 
the cloud vertical axis is 10 percent of the concentration at the cloud axis. 
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The entrainment coefficient is defined to be the ratio of cloud radius to the cloud 
centroid height. This empirical coefficient is very difficult to evaluate because of the 
complex cloud geometry and because it varies so much with height. Current entrainment 
coefficients (7 = 0.64) were obtained from cloud photographs [8,9]. 

The alongwind and crosswind source dimensions required for input to the NASA/ 
MSFC Multilayer Diffusion Models are defined for each layer by 


“xotK} =<-yo<K} 


7z/2.15 ; z < Zj^ 

7 ( 2 Zj^ - z')/2.15 > 93 m; z > Zj„ 


(5.9) 


where 


z' = midpoint of the Kth layer 
= (zbk + ztK^/2 


The quantities zjj^ and zgj^ are, respectively, the height of the top and base of the Kth 
layer. 

The corresponding vertical source dimension for each layer was calculated from 
the expression 

‘^zo "" ^^TK " 


Equation (5.10) applies to a rectangular material distribution which has been assumed to 
apply along the vertical in the Kth layer. 


E. Calculation of the Vertical Source Strength Distribution 


The fraction of material by weight in each of the K layers F{K} for the launches 
was calculated from the expression 


F{K) 


QP {^tkI ; k - 1 

Q(P {zjk} - P {zrk} ); > 1 


(5.11) 


P{zbk} is the integral of the Gaussian probability function between minus infinity and 
the base of the Kth layer zgj^, and is equal to P{zgj^ - Sigma (a) is equal to 

7{z = Zn,i)/2.15. 
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The NASA/MSFC Multilayer Diffusion Models described here require that 
source strength in each of the K layers be specified per unit height. Since the desired 
concentration units for HCI and CO are parts per million, the complete expression for the 
source strength model input for the Kth layer is 



For AI2O3, the desired concentration units are milligrams per cubic meter and the com- 
plete expression for source strength in the Kth layer is 


Qk = 


F{K)__ 
(ztk - zbk) 


10 ^ mg 


(5.13) 


Equations (5.11) and (5.12) were used to obtain the model input values of for the 
various static test results. 

F. Composition of Rocket Exhaust Effluents 

The composition of the rocket exhaust effluents varies between vehicles in accord 
with vehicle sizes and motor types. The two major rocket classes are the ones that use 
liquid and sohd rocket propellants. 

When calculating downwind concentrations from fractional weights of materials in 
the exhaust, definite uncertainties evolve. The actual amounts of elements or compounds 
available after the exhaust material combines with the ambient atmosphere is difficult to 
obtain. Factors that may cause the fractional amounts of effluents to change in the 
ground cloud are: (1) the exhaust flame evaporates thousands of gallons of deluge water 
within the flame trench and other water being sprayed on the launcher towers, (2) sig- 
nificant amounts of materials are ablated such as concrete, steel, and paint, and (3) other 
matter such as dust, ocean salt, grease, etc., are vaporized and are contained within the 
ground cloud. Subsequently, a great deal of research must be accomplished before accu- 
rate source inventory data can be made available. This is especially needed for exhaust 
ground cloud chemistry as such clouds will be composed of both solid and liquid exhaust 
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by-products. Exhaust chemistry, especially after reaction with the air and extraneous 
material, is essential for identifying initial ground cloud source composition to make 
atmospheric diffusion computations. 

G. Meteorological Inputs 

Reliable atmospheric, thermodynamic, and kinematic profiles are required to com- 
pute diffusion estimates. Consideration must be given to such factors as: (1) local clima- 
tology, (2) large scale atmospheric conditions, (3) local atmospheric conditions, (4) topo- 
graphical features, (5) land-water interfaces, and (6) exhaust source chemistry which may 
form new compounds with the ambient air, etc. One of the most important factors is 
having a sound definition of the Earth’s planetary boundary layer phenomena. This is the 
main atmospheric layer of concern when determining downwind dispersion of exhaust 
effluents. 

The meteorological inputs for diffusion modeling are as follows: (1) wind direc- 
tion, (2) wind speed, (3) standard deviations of vertical and horizontal wind, (4) humid- 
ity, (5) atmospheric pressure, (6) temperature profile data, (7) height of stable layers, and 
(8) air density [10,1 1] . Precipitation, cloud heights and types, pressure gradient condi- 
tions, and other features of the synoptic state must also be considered. The height of the 
stable layer is needed because it dictates the height the hot buoyant exhaust clouds will 
stabilize, especially when dealing with the larger vehicle exhausts. 
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CHAPTER VI 


EVALUATION OF ELECTRETS 

To evaluate the electrets, the total chlorine counts obtained from the SEM from 
the surface of the electret after exposure to pollution is a measure of the quantity of 
material collected. The counts were converted to ppm and compared with HCl measuring 
equipment and computed values from the NASA/MSFC Multilayer Diffusion Model. 

A total of 18 static test firings were conducted at MSEC Atmospheric Test Facil- 
ity 1 16. For each of the test firings, meteorological data were obtained (Tables 2 through 
19, pages 12 through 21 ) and the NASA/MSFC Multilayer Diffusion Model was used to 
predict the downwind concentration of HCl. A test matrix for the static tests is presented 
in Table 21. The coulometer, millipore filter, and bubblers used in these tests were sup- 
plied and operated by personnel from Arnold Engineering Development Center, 

Tullahoma, Tennessee. The calibration of the electrets was made from Test 18 when all 
the instrumentation appeared to assess the ground cloud for an exposure time of 10 s. 
Photographic data and personal observation indicated that the ground cloud passed uni- 
formly over the instruments during this interval. Hence, the calibration of the electrets 
was made from this test and is used in reporting the comparison of the electrets with the 
diffusion model and various HCl measuring devices. 

An evaluation of the electrets was done during the chamber tests at Arnold Engi- 
neering Development Center to verify the previously mentioned results. The electrets com- 
pared with the millipore filter. 

During the Viking I launch to Mars from KSC on August 20, 1975, an array of 
electrets was placed alongside standard HCl instrumentatiow, bubblers, and coulometers, 
0.25 km from the launch pad. During this initial flight, all power was turned off from 
T minus ten to T plus ten minutes at the launch complex. No direct comparison could be 
made with the HCl measuring equipment, but comparisons were made with the computed 
NASA/MSFC Multilayer Diffusion Model. These comparisons were made at various 
distances from the launch pad, 0.25 to 2.59 km on land and 6.80 km on the ocean. 
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Table 21. Test Matrix of MSFC Static Tests 


Test 

No. 

Date 

Time 

(CST) 

Place 


Measuring Devices 

NASA/ 

MSFC 

Diffusion 

Model 

Radiosonde 

Configuration 

Electrets 

Bubbler 

Geomet 

Coulometer 

Millipore 

Filter 

1 

August 16, 1974 

1442 

MSFC 

Solid Motor 


X 



X 

X 

X 

2 

August 30, 1 974 

1302 

MSFC 

Solid Motor 






X 

X 

3 

November 19, 1974 

2131 

MSFC 

Solid and LOX Engines 

X 





X 

X 

4 

November 22, 1974 

1830 

MSFC 

Solid and LOX Engines 

X 

X 



X 

X 

X 

5 

December 22, 1974 

1340 

MSFC 

Solid and LOX Engines 

X 





X 

X 

6 

January 7, 1975 

1759 

MSFC 

Solid and LOX Engines 






X 

X 

7 

January 9, 1975 

1758 

MSFC 

Solid and LOX Engines 


X 



X 

X 

X 

8 

January 1 1, 1975 

1435 

MSFC 

Solid and LOX Engines 






X 

X 

9 

January 1 5, 1975 

1550 

MSFC 

Solid and LOX Engines 






X 

X 

10 

January 16, 1975 

1640 

MSFC 

Solid and LOX Engines 


X 



X 

X 

X 

11 

January 18, 1975 

1450 

MSFC 

Solid and LOX Engines 






X 

X 

12 

January 21, 1975 

1 145 

MSFC 

Solid and LOX Engines 






X 

X 

13 

January 23, 1975 

1958 

MSFC 

Solid and LOX Engines 

X 





X 

X 

14 

January 27, 1975 

1822 

MSFC 

Solid and LOX Engines 






X 

X 

15 

February 8, 1975 

1254 

MSFC 

Solid and LOX Engines 


X 

X 

X 

X 

X 

X 

16 

February 12, 1975 

1700 

MSFC 

Solid and LOX Engines 






X 

X 

17 

February 22, 1975 

1528 

MSFC 

Solid and LOX Engines 






X 

X 

18 

March 25, 1975 

1700 

MSFC 

Solid and LOX Engines 

X 

X 

X 

X 

X 

X 

X 



A. Static Test Firings 


A summary of the data obtained from 18 static tests is presented in this chapter. 
From the computed NASA/MSFC Multilayer Diffusion Model, two figures are shown. The 
first figure is the computed HCl isopleths versus the crosswind and downwind distance. 
The second figure is the computed HCl maximum centerline concentration versus down- 
wind distance. Figure 4 is the HCl isopleths versus the crosswind and downwind distance 
for Static Test No. 1. 



Figure 4. Comparison of HCl concentration isopleths computed 
by NASA/MSFC Multilayer Model with measured values: 
Test No. 1 using Model 4 (height 2.0 m, azimuth 
bearing 43.0°). 
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The MSFC Test Facility 1 16 is shown in Figure 5. The Tomahawk missile is 
situated in the test stand ready for Static Test No. 1. Figure 6 illustrates a closeup of the 
rocket exhaust effluents being fired by the Tomahawk missile during Static Test No. 1. 

Static Test No. 1 

This was the first 6.4 percent SSV static test at MSFC Test Facility 116. The 
bubblers which were positioned in the superstructure were swept away by the force of 
the exhaust. The millipore filter which was 75 m from the flame trench measured a 
dosage of 410 ppm-s. The cloud passed by this detector in approximately 10 s, giving a 
concentration of 41.0 ppm [4]. Subsequent tests indicated there is a dramatic drop in 
concentration at the ground level when the cloud rises and becomes airborne. Photo- 
graphs indicated the ground cloud was airborne before 100 m. At 100 m downwind dis- 
tance, the MSFC Multilayer Diffusion Model indicated 1.2 ppm. Figure 7 shows 
the maximum HCl centerline concentration versus downwind distance. The photo 
sequence for the static test is shown in Figure 8. The movie sequence was taken with a 
35 mm Automax camera with a focal length of 50 mm and a frame speed of 2/s from 
the Saturn Static Test Stand (Building 4670), approximately 1400 m due west from Test 
Stand 116. At 80 and 90 s, the dark cloud area was due to a grass fire. 

Static Test No. 2 

The exhaust cloud from Static Test No. 2 was photographed for eight minutes 
after cloud stabilization, illustrating its travel downwind. 

The photosequence is shown in Figure 9 taken from the Saturn Static Test Stand. 
Figures 10 and 11 illustrate the HCl isopleths and HCl centerline downwind concentra- 
tion, respectively. 

Static Test No. 3 

The electrets were experimentally used for the first time in Test 3. From the 
experience gained from the two previous tests, it was determined that to obtain measure- 
ment it was necessary to get close and in line with the flame trench, which directs the 
exhaust away from the pad but not close enough to lose the electrets from the force of 
the rocket exhaust. The electret was placed 75 m in line with the flame trench. Also, 

during the previous tests it was observed from the block house (Building 4596) that the 
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Figure 5. MSFC Test Facility 116 (6.4 percent scaled model of Space Shuttle) 
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Figure 7. Computed maximum centerline HCl concentrations and 
measured values: Static Test No. 1 (height 2.0 m, 
azimuth bearing 43.0°). 
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Figure 8. (Continued). 
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Figure 10. Computed HCl concentration isopleths from NASA/MSFC 
Multilayer Diffusion Model: Static Test No. 2 using Model 4 
(height 2.0 m, azimuth bearing 64.2°). 
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Figure 1 1 . Computed maximum centerline HCl concentration from Model 4 
for Static Test No. 2 (height 2.0 m, azimuth bearing 64.2°). 
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cloud definitely passed over the electrets and that measurements should be obtained. It 
also was observed that the ground cloud hugged the ground for approximately 10 s and 
then became airborne. Later tests (15 and 18) indicated that the chemiluminescent HCl 
detector (Geomet) shows on the oscillograph the exhaust cloud passage for approximately 
10 s where measurements of 5.0 ppm and >40.0 ppm were obtained. Most importantly, 
after the exhaust cloud became airborne, the concentration dropped dramatically to zero. 
Although there were no standard measurements taken during Static Test No. 3, a com- 
parison with the NASA/MSFC Multilayer Diffusion Model indicated approximately 22.6 
ppm at 100 m. The electrets were positioned at approximately 75 m, and a measurement 
of 49.1 ppm was obtained. (See Test 18 for calibration data.) Figure 12 illustrates the 
HCl maximum centerline concentration, while Figure 13 shows the HCl isopleths versus 
downwind distance. 

Static Test No. 4 

During this test, standard HCl instrumentation and the electrets were used. The 
millipore filter and bubblers were positioned 150 m from the flame trench. Only one 
measurement was obtained by a millipore filter, 29 ppm-s, averaged over 10 s to a con- 
centration value of 2.9 ppm. Two bubblers measured zero. This firing occurred after dark 
and the visual observation indicated that the winds were very light, and the exhaust cloud 
lifted over the samples. Low ground level concentrations are obtained after the cloud 
becomes airborne. The electrets were positioned closer to the flame trench with the 
intention of obtaining a measurement. At 75 m, a measurement of 39.4 ppm was 
obtained in the ground cloud. Figure 14 illustrates the comparison of the measured values 
from the electrets, millipore filter, and bubbler with model. As shown in Figure 15, a 
measurement of approximately 79.0 ppm at 100 m was obtained from the NASA/MSFC 
Multilayer Diffusion Model. As indicated by Dumbauld et al. [5], measurements very 
close to the pad are uncertain for modeling purposes. 

Static Test No. 5 

A measurement of 5.7 ppm was obtained from the electrets positioned at 
75 m from the flame trench. Figure 16 illustrates the HCl isopleths versus crosswind 
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Figure 12. Comparison of HCI concentration isopleths computed 
by NASA/MSFC Multilayer Diffusion Model with measured 
values: Static Test No. 3 using Model 4 (height 2.0 m, 
azimuth bearing 333.6°). 
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Figure 14. Comparison of HCl concentration isopleths computed 
by NASA/MSFC Multilayer Diffusion Model with measured 
values: Static Test No. 4 (height 2.0 m, 
azimuth bearing 307.6°). 
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Figure 16. Comparison of HCI concentration isopleths computed 
by NASA/MSFC Multilayer Diffusion Model with measured 
values: Static Test No. 5 using Model 4 (height 2.0 m, 
azimuth bearing 10.4°). 
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distance and downwind distance obtained from the model. Figure 17 shows the HCl 
centerline concentration for downwind distance. At 100 m, the predicted value was 
approximately 1.7 ppm. 

Static Test No. 6 

Figure 1 8 shows the HCl isopleths versus crosswind and downwind distance and 
Figure 19 the HCl centerline concentration versus downwind distance obtained from the 
NASA/MSFC Multilayer Diffusion Model. 

Static Test No. 7 

This test was significant because the NASA/MSFC Multilayer Diffusion Model 
computed zero readings, and the instruments obtained zero readings. After deployment of 
the bubbler and millipore filter at 100 m, a delay occurred during the test. There was no 
chance to realign the instruments for the variable wind. Observations indicated that the 
positioning of the array of instrumentation missed the cloud completely, and zero values 
were indicated on the millipore filter and bubbler. Zero values were obtained on the HCl 
isopleths versus crosswind distance and downwind distance from the model as shown in 
Figure 20. Figure 21 shows the HCl maximum centerline concentration versus downwind 
distance obtained from the model. 

Static Test No. 8 

Figures 22 and 23 illustrate the predicted values obtained from the NASA/MSFC 
Multilayer Diffusion Model. The HCl isopleths versus downwind distance are shown in 
Figure 22, while Figure 23 illustrates the HCl maximum centerline concentration. 

Static Test No. 9 

Figures 24 and 25 illustrate the predicted measurements obtained from the NASA/ 
MSFC Multilayer Diffusion Model. The HCl isopleths versus crosswind and downwind 
distance are presented in Figure 24, while the maximum HCl centerline concentration is 
illustrated in Figure 25. 


58 




CROSSWIND DISTANCE (m) 





1000 10000 100000 
DOWNWIND DISTANCE (m) 


Figure 19. Computed maximum centerline HCl concentration 
from Model 4 for Static Test No. 6^ (height 2.0 m, 
azimuth bearing 324.4°). 
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Figure 21. Computed maximum centerline HCl concentration 
from Model 4 for Static Test No. 7 (height 2.0 m, 
azimuth bearing 316.4°). 
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Figure 22. Computed HCl concentration isopleths by NASA/MSFC 
Multilayer Diffusion Model; Static Test No. 8 using Model 4 
(height 2.0 m, azimuth bearing 75.0°). 
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Figure 23. Computed maximum centerline HCl concentration 
from Model 4 for Static Test No. 8 (height 2.0 m, 
azimuth bearing 75.0°). 
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Figure 25. Computed maximum centerline HCl concentration 
from Model 4 for Static Test No. 9 (height 2.0 m, 
azimuth bearing 75.8°). 
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Static Test No. 10 


From previous tests it v/as learned that to get good test results, it was necessary 
to have the instrumentation close to the launch pad. The cloud passes over the sensors in 
10 to 12 s and once the cloud becomes airborne, measurements are difficult to obtain. 
Hence, the millipore filters and bubblers were placed 75 m from and in line with the 
flame trench which directs the rocket exhaust. Figure 26 illustrates the comparison of the 
predicted HCl isopleths versus crosswind and downwind distance and the measured values 
of the HCl detectors. The millipore filters and bubbler measured from 3.2 to 10.0 ppm, 
which corresponded very well with the NASA/MSFC Multilayer Diffusion Model of 8.0 
ppm at 100 m downwind distance in Figure 27. 

Static Test No. 1 1 

Figures 28 and 29 illustrate the predicted values obtained from the NASA/MSFC 
Multilayer Diffusion Model. The HCl concentration isopleths versus crosswind and down- 
wind distances are shown in Figure 28, while Figure 29 illustrates the maximum HCl 
centerline concentration values. 

Static Test No. 12 

Figures 30 and 3 1 illustrate the predicted values of HCl obtained from the NASA/ 
MSFC Multilayer Diffusion Model. The HCl concentration isopleths versus crosswind and 
downwind distances are shown in Figure 30, while Figure 31 depicts the maximum 
centerline HCl concentration versus downwind distance. 

Static Test No. 13 

Figures 32 and 33 are the plots of the predicted values obtained from the NASA/ 
MSFC Multilayer Diffusion Model. The HCl concentration isopleths versus crosswind and 
downwind distances are shown in Figure 32, while Figure 33 illustrates the maximum 
HCl centerline concentration versus downwind distance. The electrets measurement of 
39.5 ppm at 100 m (position A) compared favorably with the models maximum center- 
line HCl concentration of 44.7 ppm. At position B, approximately 150 m from the flame 
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Figure 26. Comparison of HCl concentration isopleths computed 
by NASA/MSFC Multilayer Diffusion Model with measured 
values; Static Test No. 10 (height 2.0 m, 
azimuth bearing 171.4°). 
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Figure 28. Computed HCl concentration isopleths by NASA/MSFC 
Multilayer Diffusion Model: Static Test No. 1 1 using Model 4 
(height 2.0 m, azimuth bearing 114.2°). 
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Figure 30. Computed HCl concentration isopleths by NASA/MSFC 
Multilayer Diffusion Model: Static Test No. 12 using Model 4 
(height 2.0 m, azimuth bearing 328.6°). 
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Figure 31. Computed maximum centerline HCl concentration 
from Model 4 for Static Test No. 12 (height 2.0 m, 
azimuth bearing 328.6°). 
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Figure 32. Comparison of HCl concentration isopleths computed 
by NASA/MSFC Multilayer Diffusion Model with measured 
values: Static Test No. 13 using Model 4 (height 2.0 m, 
azimuth bearing 260.7°). 
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Figure 33. Computed maximum HCl concentration and measured 
values; Static Test No. 13 (height 2.0 m, 
azimuth bearing 260.7°). 
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trench and 20° off the direct line from the flame trench, a measurement of 1 1.3 ppm 
could be extrapolated between the isopleths of 8 and 16 ppm contour. 

S tatic Test No. 14 

Figures 34 and 35 illustrate the predicted values obtained from the NASA/MSFC 
Multilayer Diffusion Model. The HCl concentration isopleths versus crosswind and down- 
wind distances are shown in Figure 34, while Figure 35 illustrates the maximum HCl 
centerline concentration values. 

Static Test No. 15 

Test No. 15 was an early Saturday morning test. The wind direction was from 
220° and the wind speed was 2.1 m/s at the surface and at 158 m. At 325 m, the wind 
speed was 2.2 m/s from 210°. The potential temperature gradient was approximately 
0.0005 K/m. Observation after firing indicated the exhaust cloud was a little more 
buoyant as compared to the previous tests and started to rise before the end of the launch 
pad at 100 m. With the low wind speed, the cloud lingered for approximately 25 s over 
the array of detectors and then slowly dissipated. Figures 36 and 37 illustrate the pre- 
dicted values obtained from the NASA/MSFC Multilayer Diffusion Model. Figure 36 
shows the HCl isopleths versus crosswind and downwind distances, and measured values 
from the HCl detectors which compared favorably. Figure 37 depicts the maximum HCl 
centerline concentration versus downwind distance of the model. 

Static Test No. 16 

Figures 38 and 39 illustrate the HCl predictions obtained from the NASA/MSFC 
Multilayer Diffusion Model. Figure 38 shows the HCl concentration isopleths versus 
crosswind and downwind distance, while Figure 39 illustrates the maximum HCl centerline 
concentration versus downwind distance. 

Static Test No. 17 

Figures 40 and 41 illustrate the HCl prediction values obtained from the NASA/ 
MSFC Multilayer Diffusion Model. Figure 40 shows the HCl concentration isopleths 
versus crosswind and downwind distance, while Figure 41 illustrates the maximum HCl 
centerline concentration versus downwind distance. 
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Figure 36. Comparison of HCI concentration isopleths computed 
by NASA/MSFC Multilayer Diffusion Model with measured 
values: Static Test No. 15 (height 2.0 m, 
azimuth bearing 35.1°). 
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Figure 37. Computed maximum centerline HCl concentration 
from Model 4 and measured values: Static Test No. 15 
(height 2.0 m, azimuth 35.1°). 
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Figure 39. Computed maximum centerline HCl concentration 
from Model 4 for Static Test No. 16 (height 2.0 m, 
azimuth bearing 182.2°). 
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Figure 41. Computed maximum centerline HCl concentration 
from Model 4 for Static Test No. 17 (height 2.0 m, 
azimuth bearing 322.0°). 



Test Firing No. 18 

Shortly after ignition, the solid rocket motors broke loose from their mounts. 

One of the motors lodged in the test stand and the other landed in a wooded area. It was 
fortunate that all the instruments were close to the pad (approximately 75 m). Although 
it was not a normal static test firing because of the malfunction, it was a good test data 
event because the electrets were readily compared to other instruments. The average 
exposure time of the electrets and sensors to the ground cloud as observed from the 
block house and analysis of the photographs was 10 s. The dense exhaust cloud covered 
all the instruments simultaneously before it rose from the ground approximately 100 m 
from the test stand. The sensors were positioned prior to the test at 75 m. The sensors 
were placed in close because previous experience indicated once the exhaust cloud rises 
from the ground, measurements are hard to obtain. 

Figures 42 and 43 illustrate the HCl predictions obtained from the NASA/MSFC 
Mulitlayer Diffusion Model. 

An in-depth analysis is made of this test as the calibration of the electrets is 
determined from this test. 

Table 22 lists results obtained from Test 18 for the various HCl measuring instru- 
ments. Concentrations obtained using the chemiluminescent HCl analyzer from the U.S. 

Air Force School of Medicine, Brooks Air Force Base, Texas, and bubblers, coulometers, 
and millipore filters from Arnold Engineering Development Center, Tullahoma, Tennessee, 
were compared to concentrations obtained simultaneously from the electrets studies. 

Figure 44 illustrates the electrets and millipore filter positioned 75 m from the 
flame trench. Figure 45 shows the bubbler which was directly in line with the flame 
trench during the static test firing. 

Direct comparisons were made at positions 1 and 8 of HCl measurements made 
with a millipore filter and a bubbler (Table 22). The electrets were kept beside the pre- 
viously mentioned instruments to calibrate the pollution collected on the electrets sur- 
face. For example, at position 1, approximately 75 m from the flame trench and 85° 
from true north, the millipore filter measured 21.1 ppm. At the same angle and distance, 
the electrets measurement from the SEM and X-ray microprobe analysis gave 9682 
counts. Also, at position 8, the bubbler measured 66.2 ppm, and 26 204 counts were 
obtained from the SEM and X-ray microscope analysis. From the above measurements, a 
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Figure 42. Comparison of HCl concentration isopleths computed 
by NASA/MSFC Multilayer Diffusion Model with measured 
values: Static Test No. 18 (height 2.0 m, 
azimuth bearing 122.7°). 
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Figure 43. Computed maximum centerline HCl concentration 
from Model 4 and measured values: Static Test No. 18 
(height 2.0 m, azimuth bearing 122.7°). 




Table 22. Comparison of HCl Concentration Measurements Using Various 
Instruments to Electret Measurements During a 6.4 Percent Scale 
Model Test at MSFC Test Facility 116 (Test 18) 


Position 

(P) 

Angle (degree)/ 
Distance (m) 
(Measurement 
from True 
North) 

1 

Instrumentation 

Measurement 

(ppm) 

Measurement 
Obtained Directly 
from Electrets 
Counts 
(Calibrated 
Values)^ 

Estimated 
Values (ppm) 

1 

85/75 

Millipore Filter and 
Electret 

21. la 

9 682 a 

22.65 

2 

30/75 

Electret 

30.32 

(See Estimated 
Values) 

12 961 

30.32 

3 

75/75 

Bubbler 

31.3 

13 375b 

31.30 

4 

75/70 

Electret 

37.73 

(See Estimated 
Values) 

16 125 

37.73 

5 

50/75 

Millipore Filter 

45.8 

19 572b 

45.79 

6 

45/75 

Millipore Filter 

55.0 

23 503b 

54.99 

7 

65/75 

Chemiluminescent 

65.4 

27 948b 

65.39 

8 

60/75 

Bubbler and Electrets 

66 . 2 a 

26 204a 

61.31 


a. Direct comparison made between millipore filter or bubbler and electrets. 

b. Calibrated value (427.35 counts = 1 ppm). 






calibration of the counts from the electrets against the instruments (millipore filter and 
chemiluminescent analyzer) was obtained. This was found to be 427.35 counts equal 1 
ppm. It may be noted that to obtain consistent results, the SEM analyses of the electrets 
were done under identical conditions. Also, the exhaust cloud passage time over the test 
instruments was approximately 10 s, which was determined from photographs taken dur- 
ing the test. 

Figure 46 illustrates the comparison of the electrets measurements to HCl con- 
centrations from various instruments such as the millipore filter, bubbler, and chemilumi- 
nescent HCl analyzer from Test 18. 

In addition to the results in Figure 46, Figure 47 shows the distribution of con- 
centrations of the pollution from the test site at various angles for a constant distance. It 
appears that when a dense flow of air from the rocket exhaust passes an electret that is 
adjacent to a millipore filter, bubbler, etc., the reading is high; and when a less dense 
flow of air from the rocket exhaust passes, a lower measurement of concentration is 
obtained. Figure 48 is a plot similar to that of Figure 46 but shows the dosage (ppm-s) 
versus electret electron counts from Test 18. 

Figure 49 illustrates a summary of the measured data from the 1 8 static test 
firings at MSFC compared to the NASA/MSFC Multilayer Diffusion Model. It is impor- 
tant to note that the all the measured data from the electrets, millipore filter, bubbler, 
chemiluminescent detector, and coulometer were inside the upper and lower bounds of 
the diffusion model. These bounds were obtained by plotting the computed centerline 
HCl concentration from the 18 static tests and the upper and lower bounds of these tests 
encompassed the measured HCl values. 

B. Chamber Tests 

The chamber used for testing at Arnold Engineering Development Center is a 
0.914 m^ (3 ft^) chamber made of plexiglass. The sample of the rocket fuel is ignited by 
triggering an electrically heated nichrome wire. The nichrome wire ignites the solid fuel 
which then burns for approximately 1 s. The holder of the millipore filter and electrets is 
placed into the chamber before the firing sequence. Exposure times of 30, 60, 90, and 
120 s were determined using a digital timer. The mixing of the exhaust products results 
from the heat of the fuel and the resultant turbulent mixing. This lasts for several 
minutes, and observations indicated that it is a slow mixing process. 
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Figure 46. Comparison of HCl concentrations (ppm) from various instruments 
to electret electron counts during 6.4 percent SRB tests at MSFC (Test 18). 
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Figure 47. Measured distribution of rocket effluents (ppm) versus degrees from true 
north during 6.4 percent scale model testing at MSFC (Test 18). 
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Figure 48. Comparison of HCl dosages (ppm-s) from various instruments to electret 
electron counts during 6.4 percent SRB tests at MSFC (Test 18). 
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Figure 49. Comparison of measured values of HCl during 18 static tests at MSFC 
with the envelope of the upper and lower bounds obtained from the 
NASA/MSFC Multilayer Diffusion Model. 


The density or distribution of the exhaust products for the 1 s firing in a cham- 
ber is not as uniform as a typical exhaust cloud observed [7,8,9] in the atmosphere 
where there is a mass flow of air over the instruments (wind speeds of 2 to 15 m/s 
depending on the test day). By observation the distribution sensed by the millipore filters 
and electrets in the chamber is a very slow mixing process in comparison to the mass 
flow of exhaust clouds observed in rocket testing in the atmosphere. It appears, however, 
that the initial reading for both instruments in the chamber is the best reading. The flow 
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at this point is more dense and uniform. After a period of time, the distribution becomes 
spotty in parts of the chamber and the exhaust gases decay. This can be seen in Figure 
50 and Table 23. In Figure 50, the concentration for the millipore filter, after 60 s, two 
readings of 90 s, and 120 s drops off dramatically. In Figure 51, however, the electron 
counts generally increase with greater exposure time. 

From Table 24, the localized instantaneous value for the electret count of 7366 
and the millipore measurement of 20 ppm at 30s compares favorably when converted 
using 427.35 counts equal 1 ppm from the atmospheric tests. In this case, 17.25 ppm 
would be obtained from the electret measurement versus 20 ppm for the millipore filter. 

Figure 52 illustrates the electret being positioned in the chamber alongside the 
millipore filter. 
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Figure 50. Duration of exposure (s) versus concentration measurements of HCl 
from millipore filter during chamber tests at Arnold Engineering 
Development Center. 
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Table 23. Chamber Test Data HCl Concentration Measurements 
from Millipore Filter Versus Hlectret Flcctron Counts 
for Times of 30, 60, 90, and 120 s at Arnold 
Engineering Development Center 




Time of Exposure (s) 


Instrumentation 

30 

60 

90 

120 

Electret Electron Counts 

7366 

10 845 

10 263 

15 648 

Millipore Filter (ppm) 

20 

107 

263 

67 
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Figure 51. Duration of exposure (s) versus electret electron counts 
during chamber tests at Arnold Engineering Development Center. 
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Table 24. Viking I Test Matrix 


Electret 

Site 

Time 

Out 

COST) 

Time in 
(CST) 

Distance 

(km) 

Azimuth 

(degree) 

Counts 

Counts Less 
Calibration 

Counts Normalized 
for Time of 
Clouds Passage 
(counts/62.9) 

Normalized Count 

= ppm 

427.35 Counts/ppm 

NASA/MSFC 
Model 4 
(ppm) 

V, 

Tower (Prelaunch Calibration) 

1250 

1334 

0.25 

85.0 

4 076 

- 

- 

- 

- 

Vj 

Tower (Launch) 

1250 

1930 

0.25 

85.0 

9 453 

5 377 

85.4 

0.20 

0.15 

V3 

Tower (Launch) 

1250 

1930 

0.25 

85.0 

16 277 

12201 

193.97 

0.45 

0.15 

V4 

Tower (Launch) 

1255 

Lost (Blast) 

0.25 

85.0 

Lost 

- 

- 

- 

- 

Vs 

FBI 

1323 

2045 

2.53 

147.3 

37 688 

33 612 

534.3 

1.25 

1.09 

Vs 

fb: 

1348 

2008 

2.03 

170.0 

9 534 

5 438 

86.4 

0.20 

0.79 

V, 

FB3 

1310 

2037 

2.59 

344.7 

35 608 

31 532 

501.3 

1.17 

l.II 

Vs 

Tower (Postlaunch Calibration) 

1930 

1941 

0.25 

85.0 

7 786 

- 

- 

- 

- 

V, 

Boat (Prelaunch Calibration) 

1500 

1600 

6.80 

18.0 

3 1 10 


- 

- 

- 

V,o 

Boat (Vertical) 

1500 

1900 

6.80 

18.0 

5 286 

2 176 

34.6 

0.08 

1.10 

V,, 

Boat (Horizontal) 

1500 

1910 

6.80 

18.0 

24 992 

21 882 

347.8 

0.81 

1.10 

V ,3 

Boat (Post Calibration) 

1900 

2000 

6.80 

18.0 

5 940 

- 

- 

- 

- 





Figure 52. Comparison of electrets (single arrow) and millipore filter 
(double arrow) in a chamber at Arnold Engineering 
Development Center. 


C. Viking I Launch 


The Titan IIIC missile at Launch Complex 41 ready for liftoff on August 20, 
1975, for the Viking I mission to Mars is shown in Figure 53. 

It is apparent during the deployment of the instruments to measure the effluents 
from the exhaust cloud that the potential cloud trajectory is of major importance. 
Synoptic weather data from the National Weather Service Network at T minus eight 
hours indicated a nearly equal potential for both land and sea trajectory for the cloud. 
The repositioning and reevaluation of the cloud path took place at T minus four hours. 
The land cloud trajectory became more probable than the sea cloud trajectory. Haze as 
well as cloud expansion resulted in loss of visible cloud tracking after 26 min. Observa- 
tion of the cloud trajectory appeared to be 320° to 340° from the pad. 

Since the Viking Project Team cut the power off from T minus 10 minutes to T 
plus ten minutes during the Viking I launch, a direct comparison of the instruments at 
the sites could not be made with the electrets. 

The counts from the SEM as presented in Table 24 were adjusted by subtracting 
the background counts from the measurements. From the calibration of the 6.4 percent 
rocket motor firing, the exhaust cloud took an average of 10 s to travel over the sensors, 
and the calibration of 427.35 counts/ppm was based on this 10 s value. During the 
Viking I mission, the cloud took 629 s to pass over the sensor or 62.9 times more as 
compared to the 6.4 percent model on the Tomahawk missile at MSFC. A comparison of 
the results of the electrets and the NASA/MSFC Multilayer Diffusion Model is given in 
Table 24. 
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Figure 53. Titan IIIC poised at Launch Complex 41 KSC for launching 
of Viking I mission to Mars on August 20, 1975. 
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CHAPTER VII 


CONCLUSIONS 

The purpose of this research investigation was to evaluate the electrets, a new 
detecting contamination device used to assess the chemical compositions of rocket efflu- 
ents. In assessing the effectiveness of the electrets, the following was concluded: 

The electrets have been used successfully in collecting rocket exhaust effluents. It 
was possible by the use of the SEM and X-ray microprobe analysis to obtain spectra, 
giving the electret a multipollutant capability. In addition, simplicity of deployment and 
speed of assessment and analysis make the electrets a valuable tool in measuring pollu- 
tants from a space vehicle rocket exhaust. 

The electret compared favorably with other HCl measuring devices. It is a valuable 
complementary device to other measuring systems. Summary of the measured data from 
the electrets and other HCl detectors is within the upper and lower bounds of the com- 
puted HCl concentrations using the NASA/MSFC Multilayer Diffusion Model. 

Future research should be designed for extensive investigations to obtain quantita- 
tive results. Ideally, the measurements should be repeated so that reliability of the results 
could be assured. Continued chamber and field tests should achieve these results. 
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APPENDIX A 


CALCULATION OF THERMODYNAMIC VARIABLES 
FROM RAWINSONDE DATA 


The equations used for calculation of thermodynamic variables from measure- 
ments of altitude, temperature, and relative humidity obtained from the GMD-4, AMQ-9 
rawinsonde system are summarized herein; these equations, originally developed for the 
GMD-2 system [11], must be used in conjunction with the list of symbols and units pro- 
vided at the end of this appendix. 

Atmospheric Density, p 
p 

p = 348.38 — (A-1) 

T 

*v 


Pressure, P 

P = P' 10 -(h-h')/( 221 . 266 Tyj„) 

Geopotential Height, h 

^o ^e^ 

h = 

9.8 rg + H 

Virtual Temperature, Ty 

Ty = T(l +0.376932 e/P') 

Mean Virtual Temperature, T^j^^ 

T' +T 
V ‘v 

X = 

vm 0 


(A-2) 


(A-3) 


(A-4) 


(A-5) 


Vapor Pressure, e 

e = 6.11 fj) io7-5t/(t+237.3) 


(A-6) 


Dew Point Temperature, t^j 

_ 237.3 log e- 186.527 
^36 - log e 


(A-7) 
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Potential Temperature, 0 


1000 “-2 8 8 

0 = T 

P 


Virtual Potential Temperature, 0y 
10000-2 8 8 

Absolute Humidity, 

p^ = 216.7 e/P 


Microwave Refractive Index, n 


n= 1 + 


77.6P- lle + 


374808e 


T 




10 


r6 


For data tabulation, use; 
N = (n- DIO* 

Speed of Sound, Vg 


(A-8) 


(A-9) 


(A- 10) 


(A-11) 


(A- 12) 



LIST OF SYMBOLS 


L 

/ 


Symbol 

e 

So 

h 

h' 

H 


H' 


n 

N 

P 


P' 


t 

T 

td 

Ty 

n 


vm 


’ s 
P 

Pw 

d 


Definition 

Vapor pressure (mb) 

Relative humidity expressed as a decimal 

Acceleration of gravity at geographical location of the rawinsonde 
station (m/s^ ) 

Geopotential height at the top of the layer bounded by h and h' 

(ft) 

Geopotential height at the bottom of the layer bounded by h and h' 
(ft) 

Geometric altitude at the top of the layer bounded by H and H' 

(ft) 

Geometric altitude at the bottom of the layer bounded by H and 
H' (ft) 

Microwave refractive index 

Unit of refractive index used for simplification of data tabulation 

Pressure at geopotential height h (mb) 

Pressure at geopotential height h' (mb) 

Radius of the Earth (ft) 

Temperature (°C) 

Temperature (K) 

Dew point temperature (°C) 

Virtual temperature at geopotential height h (K) 

Virtual temperature at geopotential height h' (K) 

The mean virtual temperature of layer bounded by h and h' 

(K) 

Speed of sound (kn) 

Atmospheric density (gm/m^ ) 

Absolute humidity (gm/m^) 

Potential temperature (K) 

Virtual potential temperature (K) 
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